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Fermi level tuning and a large activation gap achieved in the topological insulator 

Bi 2 Te 2 Se by Sn doping 
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We report the effect of Sn doping on the transport properties of the topological insulator Bi2Te2Se 
studied in a series of Bi2- a; Sn a; Te2Se crystals with < x < 0.02. The undoped stoichiometric 
compound (a; = 0) shows an n-type metallic behavior with its Fermi level pinned to the conduction 
band. In the doped compound, it is found that Sn acts as an acceptor and leads to a downshift 
of the Fermi level. For x > 0.004, the Fermi level is lowered into the bulk forbidden gap and the 
crystals present a resistivity considerably larger than 1 ficm at low temperatures. In those crystals, 
the high-temperature transport properties are essentially governed by thermally-activated carriers 
whose activation energy is 95-125 meV, which probably signifies the formation of a Sn-related 
impurity band. In addition, the surface conductance directly obtained from the Shubnikov-de Haas 
oscillations indicates that a surface-dominated transport can be achieved in samples with several 
/^m thickness. 

PACS numbers: 73.25.-t-i, 74.62.Dh, 72.20.My, 73.20.At 
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I. INTRODUCTION 

In a three-dimensional (3D) topological insulator (TI), 
a band inversion due to strong spin-orbit coupling in- 
duces gapless surface states (SS) consisting of spin helical 
Dirac FermionSji~— which are expected to give rise to a 
number of topological quantum phenomena^ 7 - Whereas 
such surface sensitive probes as angle-resolved photoc- 
mission spectroscopy (ARPES) and scanning tunneling 
microscope have been successfully applied to investigate 
the topological SS, 8 ~ 19 the transport study of the SS re- 
mains a challenge due to the presence of parallel bulk 
conducting channel that usually dominates the trans- 
port properties i^2r— Recently, the ternary tetradymite 
Bi2Te2Se (BTS) has become a prototype TI material for 
studying the peculiar spin and charge transport of the 
SS ) 29 i 30 because carefully prepared BTS crystals not only 
present a very low bulk conduction, but also show clear 
Shubnikov-de Haas (SdH) oscillations that reflect a high 
mobility of surface Dirac electrons. Given the lack of a 
truly bulk-insulating state in any of the known TI ma- 
terials so far,— ~— further improvement of BTS is worth 
pursuing. 

The first synthesis of BTS was made while optimiz- 
ing the continuous solid solutions between Bi2Te3 and 
Bi2Se3 for thermoelectric applications^ Nevertheless, 
unlike other members in the Bi2Te3_ a; Se 3 ; (0<x<3) fam- 
ily, the Te and Se atoms in BTS occupy distinct crystal- 
lographic sites, forming quintuple layers arranged in the 
sequential order Te-Bi-Se-Bi-Te along the c axis. 35 Such 
chalcogen ordering is believed to provide structural ba- 
sis for reducing bulk carriers and achieving high surface 
mobility 29 Unfortunately, similar to Bi2Se3, the BTS 
crystals grown from the stoichiometric melts show n-type 
metallic conduction with the electron density n e of ~10 19 
cm~ 3 . 36 When the chalcogen stoichiometry is altered in 
the starting composition, one obtain a 29 ' 36 metallic or in- 
sulating crystals at different positions along the boule, 



which is due to the inevitable phase separation during the 
solidification process, according to a newly established 
phase diagram of the Bi 2 Te3-Bi 2 Se3 system^ 7 - Therefore, 
it is desirable to explore an alternative route to obtain 
BTS crystals with a large bulk resistivity while keeping 
the chalcogen stoichiometry, which would be useful for 
achieving a high surface-carrier mobility. 

In this paper, we show that the above objective can 
be accomplished by hole doping through a Bi-site sub- 
stitution, which has already been proven to be effective 
in tuning the carrier type and density for both of the 
binary end members, Bi 2 Te 3 and Bi 2 Se3i 10 ' 22 ' 23 i 28 It is 
found that, in analogy to its role in Bi2Te3 ) 10 ' 38 ~ — the 
group IV element Sn acts as an acceptor in BTS, which 
allows us to tune the Fermi level (E-p) of this material. 
In particular, the Big-xSna; Te2Se crystals with x > 0.004 
show low-temperature resistivity reaching several f^cm, 
indicating that E-p is tuned into the bulk band gap; their 
transport properties at high temperatures signify a large 
activation gap, which probably comes from a Sn-related 
impurity band. In about a half of those samples, we 
observed SdH oscillations originating from the topolog- 
ical surface state below 14 T. Our analysis of the SdH 
oscillations gives direct evidence that one can achieve 
a surface-dominated transport in a bulk Sn-doped BTS 
crystal with a thickness of several /im. 



II. EXPERIMENTAL DETAILS 

The single crystals of Bi2_ ;E Sn a: Te2Se were grown by 
melting high purity elemental shots of Bi (99.9999%), 
Sn (99.99%), Te (99.9999%), and Se (99.999%) with a 
nominal ratio of Bi:Sn:Te:Se = (2-x):a;:2:l (a; = 0, 0.002, 
0.004, 0.006, 0.01, 0.02) in sealed evacuated quartz tubes 
at 850 °C for 48 h with periodically shaking to ensure ho- 
mogeneity, followed by cooling slowly to 500 °C and then 
annealing at that temperature for 4 days. The resulting 



crystals are easily cleaved along the (111) plane, revealing 
a shiny mirrorlike surface. The x-ray diffraction (XRD) 
analysis, which was performed on powders obtained by 
crushing the crystals, confirmed the samples to be single 
phase with chalcogen-ordered tetradymite structure. 

For transport characterizations, the crystals were cut 
into bar-shaped samples with the typical thickness of 100 
/jm after they were checked to be single domain by x-ray 
Laue analysis. The electrical leads were attached to the 
samples using room-temperature-cured silver paste in a 
six-probe configuration. The in-plane resistivity p xx and 
the Hall coefficient i?n were measured in a Quantum De- 
sign Physical Properties Measurement System (PPMS-9) 
down to 1.8 K. The magnetic field was applied along the 
C3 axis which is perpendicular to the cleaved surface. 
For each x value, the data were taken on several crys- 
tals obtained from different parts of the boule in order 
to check for sample-to-sample variation. In addition, se- 
lected samples were cleaved down to a few micrometers 
by using Kapton tapes and then brought to a 14-T mag- 
net for the detection of SdH oscillations using an ac mea- 
surement technique, in which two lock-in amplifiers were 
employed to collect the signals in both the longitudinal 
(Pxx) and transverse (p yx ) channels simultaneously. The 
measurements were carried out by sweeping the magnetic 
field between ±14 T at the rate of 0.3 T/min, during 
which the temperature was stabilized to within ±5 mK. 



III. RESULTS 
A. Resistivity 

Figure 1 (a) shows typical data for the temperature de- 
pendences of p xx in Bi2_ a; Sn a ;Te2Se crystals with differ- 
ent x values, together with the data for an insulating sam- 
ple grown from the starting composition Bi2Te1.g5Se1.05 
which we abbreviate as BTSi.os- The undoped stoichio- 
metric BTS crystal (x = 0) shows a metallic behavior in 
the whole temperature range; correspondingly, the Hall 
coefficient Ru [inset of Fig. 2(a)] is negative and nearly 
temperature independent, giving an electron carrier den- 
sity of ~1.5xl0 19 cm~ 3 at 1.8 K. This result is in agree- 
ment with the previous study^ which indicated that E-p 
is pinned to the conduction band for this composition. 
Upon doping with Sn, a drastic change in both the re- 
sistivity value and its temperature dependence was ob- 
served. It is worth noting that an insulating behavior 
of p xx is already established even for 0.1% of Sn doping 
to the Bi site (x = 0.002). For x > 0.004, the p xx val- 
ues reach a few 17cm, which are usually larger than that 
achieved in BTSi.05^ 

It should be noted that there is a marked difference 
in the p xx (T) behavior between the Sn-doped BTS (x 
> 0.004) and the BTS1.05 samples. Firstly, in the high- 
temperature region, p xx in Sn-doped BTS samples in- 
creases more steeply with decreasing temperature, as can 
be seen more clearly in the Arrhenius plot shown in Fig. 
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FIG. 1: (Color online) (a) Temperature dependences of p X x 
for Bi2- a; Sn a ;Te2Se single-crystal samples with various x. 
Note that the data with x = and x — 0.002 have been 
magnified by a factor of 500 and 10, respectively. For com- 
parison, the data for an insulating sample grown from the 
starting composition Bi2Te1.95Se1.05 is also included (dashed 
line), (b) Arrhenius plot of the p X x(T) data for the temper- 
ature range between 100 and 300 K. (c) p xx (T) data for four 
representative samples with x = 0.006 obtained from differ- 
ent parts of a boule. The data are labeled in numerical order, 
where SI and S4 represent first-to-freeze and the end parts, 
respectively. Note that the vertical axis is in logarithmic scale. 



1 (b) . The magnitude of p xx at room temperature in Sn- 
doped BTS samples varies from 80 to 110 mJ7cm, which 
are considerably larger than that in BTSi.05^ Also, the 
high-temperature slope in the Arrhenius plot, which is 
almost doping independent, corresponds to an activation 
gap A of approximately 120 meV; this is nearly three 
times larger than that found for the BTS1.05 sample (A ~ 
45 meV) . Secondly, in the low temperature region, while 
Pxx{T) of the BTS1.05 sample becomes nearly flat, the Sn- 
doped BTS samples show a more strongly temperature- 
dependent p xx (T) behavior with a weak upturn below 
~10 K, whose origin is not clear at present. 

In addition to the above observations, the Sn-doped 
BTS samples were found to show a higher degree of ho- 
mogeneity. As an example, Fig. 1(c) shows the p xx (T) 
data for four samples obtained from different parts of 
a boule for x — 0.006. One can see that all the sam- 
ples show qualitatively similar temperature dependence 
of p xx . Although the p xx values still vary between sam- 
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FIG. 2: (Color online) (a) Temperature dependences of the 
low-field Hall coefficient Ru for a series of Bi2- :E Sn ;E Te2Se 
samples, together with the data for the BTS1.05 sample. The 
inset shows a magnified view of the data with x — and x — 
0.002. (b) Arrhenius plot of Rn(T) at high temperature for a 
series of samples except for x = and 0.002 where no activa- 
tion behavior was observed. Dashed lines are the Arrhenius- 
law fittings to extract the activation energy A*, (c) Ru{T) 
data for two different samples with x — 0.004, showing both 
n-type and p-type behavior, (d) The density of introduced 
holes (Arih — n(x) - n(0)) by Sn doping plotted as a function 
of the density of Sn dopant atoms. The estimation of the car- 
rier density is based on a one-band model for simplicity. The 
dashed line corresponds to the ideal situation where each Sn 
atom donates one hole. 



pies by a factor of 4, it is notable that all of them are 
larger than 1 £7cm at low temperature. 



B. Hall coefficient 

To further investigate the doping effect of Sn in BTS, 
we measured the Hall resistivity p yx . Figure 2(a) shows 
the temperature dependences of the Hall coefficient Rn 
for the same series of BTS1.05 an d Bi2- x Sn a; Te2Se sam- 
ples; here we follow our previous definition of Ru as Rr 
= Py X /B near B = 0£L It turns out that 0.1% of Sn dop- 
ing results in an increase in |.Rh| by nearly one order of 



cm , implying that one Sn atom substituted for a Bi 
atom in BTS introduces approximately one hole, which 
is much the same as that in the case of n-type Bi2Te3. 10 
While i?H is negative for x = and 0.002 [Fig. 2(a) in- 
set], it becomes positive at x = 0.004 [Fig. 2(a) main 
panel], suggesting that an n-to-p type transition occurs 
as a result of Sn doping between x = 0.002 and 0.004. 
In fact, in the case of x = 0.004, both n- and p-type 
samples with a large |-Rh| were found in the same batch 
[Fig. 2(c)], placing this composition at the verge of such 
a transition. Somewhat unexpectedly, Rn becomes neg- 
ative again at a higher doping of x = 0.02 [Fig. 2(a)], 
whose origin is not clear. In the following, we focus on 
the temperature dependence of Ru for the samples with 
x > 0.004. 

Above ~200 K, i?H of those samples shows a thermally 
activated behavior, indicating that E-p is lowered into the 
bulk band gap. The effective activation gap A*, which 
is obtained from the Arrhenius plot of Rn(T) shown in 
Fig. 2(b), is 115, 95, 125, and 110 meV for the samples 
with x = 0.004, 0.006, 0.01, 0.02, respectively. Actually, 
the A* values for all the measured crystals fall within 
95-125 meV for those compositions. These A* values 
are not far from the A values derived from the p xx {T) 
data, and their difference is likely due to the temperature 
dependence of the carrier mobility. More importantly, 
the A* values are much larger than that found in the 
BTS1.05 sample (~65 meV), implying that the Sn doping 
brings Ep closer to the middle of the band gap. At low 
temperatures, the \Rn\ values of those Sn-doped samples 
become large, sometimes exceeding 1000 cm 3 /C which 
would correspond to the carrier density of only 6 x 10 lD 
cm -3 in a one-band model. 



C. Surface quantum oscillations 

Although Sn doping is expected to introduce impurity 
scattering, it turned out that the surface mobility re- 
mains reasonably high so that SdH oscillations can still 
be observed. As a matter of fact, traces of SdH oscilla- 
tions were detected in nearly 50% of the Sn-doped BTS 
samples measured in the 14-T magnet. Among those suc- 
cessful cases, the data taken on an n-type sample with 
x = 0.004 (thickness t — 6 nm)^ showed the simplest 
pattern of the oscillations. In the following, we present 
the analysis of this simplest case. 

Figure 3(a) shows the magnetic-field dependence of 
the transverse conductivity G yx of this sample at 1.6 
K, which was calculated from p xx and p yx . SdH oscilla- 
tions are already visible in the raw data for magnetic field 
above 10 T. After removing a smooth background, one 
can clearly see that the oscillatory part of G yx , AG yx , 
exhibits periodic maxima and minima as a function of 
1/B [Fig. 3(a) inset], establishing the existence of a well 
defined Fermi surface (FS). In Fig. 3(b), we plot the 
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FIG. 3: (Color online) (a) Transverse conductance G yx for an 
n-type sample with x = 0.004 at 1.6 K plotted as a function of 
magnetic field B applied along the C3 axis. The inset shows 
the oscillatory component of G yx , AG SI , plotted as a function 
of 1/B. The dashed lines are the guide to the eyes, (b) 
Landau- level fan diagram for oscillations in AG yx - Maxima 
and minima correspond to n + 1/4 and n + 3/4, respectively. 
The linear fitting with the slope of F — 116 T fixed by the 
oscillation frequency intersects the axis at /3 = 0.4 ± 0.1. 
(c) Dingle plot for the oscillations in AR XX , yielding Td — 
12.5 K; inset shows the temperature dependence of the SdH 
amplitudes. The solid line represents the fitting by the LK 
theory, giving m c = 0.13m e . 



1/B values corresponding to the maxima (closed circles) 
and the minima (open circles) of AG yx as a function of 
the Landau level index n, following the index assignment 
scheme in Ref. 42. From the linear fitting of the data 
with the slope fixed at the oscillation frequency obtained 
from the Fourier transform (F = 116 T), we obtain a 
finite intercept j3 — 0.4 ±0.1. Since the slope is fixed in 
this analysis, the error in /3 is relatively small; the main 
source of the error is the uncertainty in determining the 
positions of maxima and minima in the data shown in 
the inset of Fig. 3(a), and the error of ±0.1 in j3 is a 
conservative estimate. The obtained j3 of 0.4 ± 0.1 is 
reasonably close to the value j3 = 0.5 expected for mass- 
less Dirac Fermions, which points to the topological SS 
origin of the SdH oscillations. Using the Onsager relation 
F = (hc/2ire)A, where A is the extremal FS cross-section 
area, we find the Fermi wave vector fcp = 5.9 x 10 6 cm -1 , 
which corresponds to the surface carrier density n s = 2.8 
x 10 12 cm -2 for a spin-nondegenerate surface state. It 
is worth noting that if one assumes that the SdH oscilla- 
tions originate from a bulk FS, the bulk carrier density 



implied by F is of the order of 10 18 cm -3 , which is to- 
tally inconsistent with the large R~h value observed at 
1.8 K (~ -100 cm 3 /C). Hence, one can conclude that the 
SdH oscillations are certainly coming from the surface. 
From the obtained kp value and the Fermi velocity v-p ~ 
4.6 x 10 5 m/s^ Ey is estimated to be ^170 meV above 
the Dirac point, pointing to the electron character of the 
surface carriers. 



By fitting the temperature dependence of the oscilla- 
tion amplitude with the standard Lifshitz-Kosevich (LK) 
theory [inset of Fig. 3(c)], we obtained the cyclotron 
mass m c — 0.13m e , where m e is the free electron mass. 
Once m c is known, the Dingle plot [Fig. 3(c)] gives the 
Dingle temperature Td = 12.5 K, which corresponds to 
the surface quantum mobility /if dH of ^1300 cm 2 /Vs. 
This fig dH value is roughly twice as large as that reported 
for thick BTS1.05 samplesi 2 ^ Surprisingly, according to 
these results, the ratio of the estimated surface conduc- 
tance to the total conductance, given by G$/crb, is cal- 
culated to be about 1.8 in this sample; here, the surface 
sheet conductance G, = en,/^ dH ~ 5.8 x 10~ 4 fi _1 , the 
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measured overall conductivity a was 0.64 f2 crrT 
1.6 K, and t = 6 /1m. 4 - Obviously, it is unphysical to 
have G s /at > 1, and this result implies that the actual 
surface transport is hindered by steps between terraces 
created upon cleaving, giving the effective surface con- 
ductance that is smaller than the intrinsic one given by 
n s and /j, s . Note that a similar problem was reported in 
Ref. [i^. In any case, the SdH oscillations give direct esti- 
mate of the surface transport parameters, which indicate 
that a surface-dominated transport can be achieved in 
Sn-doped BTS crystals with several /mi thickness. The 
remaining bulk contribution is most likely the degenerate 
transport through the impurity band which is located in 
the bulk band gap; to further reduce the bulk contribu- 
tion, one needs to achieve the Anderson localization in 
this impurity-band transport. 



Apart from this simplest case, the patterns of the ob- 
served SdH oscillations were complicated. One of such 
examples, which was observed in a sample with x — 0.01, 
is shown in Fig. 4. As can be seen in the lower inset, 
the Fourier transform of AG yx (B~ 1 ) reveals two well- 
resolved peaks at F\ = 70 T and F^ = 132 T, reflecting 
the beating visible in the oscillation data (upper inset). 
The fcp values calculated from the Onsager relation are 
4.6 x 10 6 and 6.3 x 10 6 cm -1 for F\ and F 2 , respectively. 
These correspond to n s of 1.7 x 10 12 and 3.2 x 10 12 
cm 1 respectively. Unfortunately, as was pointed out in 
Ref. l28J, the multi-component nature of the oscillations 
prevents us from reliably extracting the cyclotron mass 
or the Dingle temperature for each component. Also, the 
Landau-level fan diagram is not very reliable for extract- 
ing j3 in the multi-component case. 
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FIG. 4: (Color online) G yx for a sample with x = 0.01 at 1.6 
K plotted as a function of B applied along the C3 axis. The 
upper inset shows the oscillatory component AG yx plotted as 
a function of 1/B. The lower inset shows the Four transform 
of AG yx (B~ ), revealing two prominent frequencies. 



Another, more plausible possibility is that, at high Sn- 
doping concentration, the wavefunctions of the Sn accep- 
tors overlap significantly, leading to the formation of an 
impurity band (IB)^^S Within this picture, increasing 
Sn content results in an increase in the density of states 
(DOS) of the IB instead of hole doping. In Sn-doped 
Bi2Te3, this Sn-related IB was shown to be located at 15 
meV below the top of the upper valence bandit How- 
ever, in the case of Sn-doped BTS, this Sn-related IB is 
most likely located within the bulk band gap, because 
the activation energy was found to be much larger in Sn- 
doped BTS than in BTS1.05, which naturally points to 
the appearance of a new IB to pin the chemical poten- 
tial in Sn-doped BTS. Note that, while there may also 
be IBs due to Se vacancies and Bi/Te antisite defects in 
Sn-doped BTS, it is most likely that E-p is pinned to the 
Sn-related IB due to its large DOS. The fact that the 
activation energy of ^120 meV is essentially unchanged 
for a range of Sn concentrations (x — 0.004-0.02) is also 
consistent with this picture. Further studies are called 
for to clarify the details of the IBs in Sn-doped BTS. 



IV. DISCUSSIONS 



CONCLUSION 



From the above results, it is clear that Sn acts as an 
acceptor in BTS. Ideally, one would expect that the total 
number of holes introduced by Sn doping, Ann , increases 
linearly with increasing x. However, as can be seen in Fig. 
2(d), Ant becomes nearly constant for x > 0.004, which 
deviates significantly from the ideal situation. Therefore, 
there must be some additional effects that lead to deac- 
tivation of Sn dopants. In this respect, a similar problem 
was noted in previous studies of Sn doping in T5i 2 Te3r^ r — 
and two possibilities have been proposed to explain the 
apparent discrepancy. 

One possibility considered for Bi2Te3 is that a part 
of the Sn atoms are built into the lattice in such a way 
that a seven-layer lamellar structure, i.e. Te-Bi-Te-Sn- 
Te-Bi-Te, is formed^ If a similar structure exists in Sn- 
doped BTS as well, it would be Se-Bi-Te-Sn-Te-Bi-Se. It 
is important to note that, in contrast to those occupy- 
ing the Bi site, the Sn atoms in this structure do not 
bring any charge to the lattice, and thus the discrepancy 
is reconciled. However, no additional diffraction peaks 
except for those corresponding to the chalcogen-ordered 
tetradymite structure were observed in the XRD data 
even for the sample with the highest Sn-doping concen- 
tration (data not shown), suggesting that the seven-layer 
lamellae, if exist, are randomly distributed in the lattice. 



We performed a systematic study of the transport 
properties of a series of Bi2_ a; Sn a ;Te2Se single crystals 
with < x < 0.02. It is found that Sn behaves as an ac- 
ceptor, which enables us to tune the Fermi level that is lo- 
cated in the conduction band in the undoped stoichiomet- 
ric compound. For > 0.004, Ep is successfully tuned into 
the bulk band gap, and the resistivity becomes as large 
as several S7cm at low temperatures. The transport prop- 
erties at high temperatures show a thermally activated 
behavior with a large activation gap, which is probably 
related to the formation of a Sn-related impurity band. 
The analysis of the SdH oscillations observed in a 6-/xm 
thick sample indicates that a surface-dominated trans- 
port can be achieved in Sn-doped BTS single crystals 
with several /jm thickness. This, along with the large ac- 
tivation gap, makes the Sn-doped BTS system well suited 
for future applications of topological insulators. 
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